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Abstract 

We discuss the group theory relevant to the ground-state baryons in large N c QCD. 
For very large representation, the group generators become classical variables. We 
find the form of the classical generators for the completely symmetric N index repre- 
sentation of SU(m) as N — > oo and derive an integral formula for the matrix elements 
of an arbitrary polynomial in the group generators between low-spin baryon states in 
the large N limit. 
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The idea of replacing the SU(3) gauge symmetry of QCD with an SU (N) symmetry and studying 
the A" — > oo limit, as articulated by 't Hooft in [|IJ, has led to an important qualitative understanding 
of some of the properties of QCD, such as Zweig's rule and the narrowness of resonances. Witten 
[T2J later provided a conceptual framework to include baryons in the theory. In Witten's picture, 
a baryon is described by a Hartree-Fock equation with each quark moving in the mean potential 
generated by all the other quarks in the baryon. Using this description for the baryons, [|3], [|. |3j 
(see also 0) one can show that to leading order in N, the low-spin baryons have a spin- flavor 
symmetry which we will denote as SU(m), e.g. for an A/-flavor theory m = 2Nf. In this letter, 
we study the group theory associated with the large representations of this SU{m). In doing so, 
we develop an elegant integral formula for matrix elements of SU(m) generators between low-spin 
baryon states. This formula provides an insight into the nature of the large A" enhancement of these 
matrix elements. 

An A" quark baryon lives in the completely antisymmetric representation of the color SU(N). 
An s-wave ground state, to satisfy fermi statistics, must then be completely symmetric under the 
SU(m) symmetry. Therefore, the representation of SU(m) relevant to large A" baryons is the 
completely symmetric combination of A" defining representations. We will examine the matrix 
elements of the group generators in this representation, which we will denote by Tg , to leading 
order in N, and will discover that many properties of the generators can be obtained rather simply 
in this limit. 

Our primary tool will be the fact that for large N, the commutator of two group generators is 
lower order in A" than the product. The product is order A^ 2 , while the commutator is order A" 
(times the structure constants, which depend on m but not on N). Thus to leading order in N, the 
generators can be simultaneously diagonalized — they become essentially classical variables @ 0- 

Explicitly calculating the traces in the representation space of products of generators and taking 
the N —>■ oo limit leads to the following form for the generators:^ 

T$\u)=T?(u) \u) 

(1) 

1 We will present a detailed proof for this result in a future article. The same result also emerges in the geometric 
quantization of SU(m). || 
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where u is an m component complex unit vector and up = u^*. Notice that this form for the 
generators has a single non-degenerate eigenvalue ((m — l)/m) and (m — 1) degenerate eigenalues 
(— 1 for the eigenvectors orthogonal to u) and that moreover the overall phase of u does not matter. 
The generators in the large N limit belong to the classical coset space, 

SU(m)/[SU(m — 1) x U(l)] = CP m_1 , (2) 

but (H) is stronger than (0), because the normalization is also fixed. For an SU(2) symmetry, this 
limit reproduces the familiar statement that angular momentum becomes a classical vector when 
the spin grows very large, since then SU(2)/U(1) is a two-sphere that labels the direction of the 
classical angular momentum. 

We show in the following section how the form for the generators, (|I]), naturally arises in the 
m dimensional harmonic oscillator. The harmonic oscillator energy eigenstate live in precisely the 
same representation of SU(m) as do the ground state baryons in large iV QCD. In section |2|, we 
find the eigenstates of the group generators, and in section [| we use these eigenstates to find an 
integral formula for the leading contribution to matrix elements of an arbitrary polynomial in the 
generators between low-spin states. 

1 m Dimensional Harmonic Oscillator 

We begin by analyzing the large N limit of the m-dimensional harmonic oscillator. In addition to 
providing a familiar system which has a well-known classical limit, this example will prove useful 
later when we attempt to construct the low-spin baryons. In the large N limit of QCD, the important 
behavior of the quarks comprising a baryon can be represented by colorless, bosonic creation and 
annihilation operators, as has been applied in the "quark representation" of and in the analysis 
of || . We will use the energy eigenstates of the harmonic oscillator as an intermediate step bridging 
between the baryon states and the SU(m) generator eigenstates. 

The m-dimensional harmonic oscillator corresponds to the Hamiltonian (setting the spring con- 
stant and the mass to 1) 

H= 1 -(^+^)=^a a + ^ (3) 
where a a ^ and ap represent the raising and lowering operators 



Besides the energy, H, the classical harmonic oscillator has two additional conserved quanities: the 
angular momentum, 

M afB = x a p f3 -x f3 p a (5) 

and the traceless symmetric tensor 

S°0 = x a + p a / - — 5 al3 H. (6) 

m 

These conserved tensors correspond, respectively to the antisymmetric and the symmetric parts of 
the SU(m) generators, 

Tp = a a ^a p - -5 Q V f a 7 = \s a ? - % - M a(S . (7) 

777 Zi Zi 

The raising operators are tensor operators transforming under SU(m) like the defining m di- 
mensional representation. The quantum states are constructed by acting on the ground state with a 
product of raising operators. Since the raising operators commute, a state with N raising operators 
transforms as the completely symmetric SU (m) tensor with N indices — the same representation in 
which arose for the large N baryons. We would expect that for large N, the conserved quantities 
would approach their classical values. Let us show that this happens. 

The most general classical motion of the m dimensional harmonic oscillator is 

x(t) = v\ cos(t) + v 2 sin(t) p(t) = —v-y sin(t) + v 2 cos(t) (8) 
The real vectors V\ and v 2 satisfy 

vl + vl = x 2 + f = 2E. (9) 

In the state with principle quantum number N, the energy is N + m/2, or approximately iV for iV 
large. Thus, in the large N limit, 

vt + v 2 2 = 2N. (10) 

The traceless tensor S a ^ becomes 

v? v{ + < v$-— 5 aP ( vl + vf] (11) 
while the angular momentum tensor becomes 

uj>£-u£uf. (12) 
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Now introduce the complex vector 

V2N u = v\ + i v 2 

so that 

u u — I 

In terms of this complex vector u, we discover that 

T% = - S a/3 — - M al3 = N (u a u p * - — 5 aP ] 
13 2 2 V m J 

which is just ([[]). 

2 The Eigenstates of the Generators 

Our ultimate goal is to describe the leading N behavior of matrix elements of the SU(m) generators 
between low spin baryon external states. The simplest states, however, with which to describe the 
action of the generators are the eigenstates introduced in equation ([[]) since on these eigenstates the 
group generators act diagonally. To connect these two pictures, we must learn how to express the 
baryons in terms of these generators. Since the set of eigenstates is complete, it should be possible 
to write any state, such as that for a baryon or for a harmonic oscillator energy eigenstate, as some 
distribution over the space labeling the eigenstates of the generators. Conceptually, it is simplest 
to proceed in two steps — first to express the harmonic oscillator energy eigenstates in terms of the 
|it)'s and then to write baryons in terms of these energy eigenstates. 

To begin, we choose a system of coordinates to describe an eigenstate of the SU(m) generators. 
Up to an overall phase, to each generator corresponds a unique unit vector, u; we shall use this 
arbitrariness of phase to set the m th component of u to be purely real. The remaining components 
of u can be expressed in terms of (m — 1) positive numbers w a and corresponding phases 9 a . Thus 
the vector u(w, 6) becomes 

u a = e tda Jw^ [a — 1, . . . , m — 1] 

(16) 

Since these coordinates are continuous, the eigenstates require a continuum normalization, 

m— 1 

(u(w',0')\u{w,6)) = J] 5{w' a -w a )8{9' a -e a ) (17) 

a=l 
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(13) 
(14) 

(15) 



which has the correct 8(w' — w) dependence to ensure that the measure is SU{m) invariant. 

Having chosen this coordinate system, we next establish the connection between the harmonic 
oscillator energy eigenstates and the u(w, 9) eigenstates. The oscillator has correctly normalized 
eigenstates of the form 



m_ r t ]n 



\n) 



where N = J2 a n a goes to infinity. Since the eigenstates are complete, we can generally write 
\n) as some distribution over the space of eigenstates \u(w, 9)), specified by a "spectral function" 
C(n, w,9): 

\n)= J dm9C(n,w,9)\u(w,9)) +OQ./N). (19) 

Acting upon both sides of this expansion with the group generators the group generators, Tp, 
determines the behavior of the function C(n, w, 9): 



(a a ^ dp — — <52a 7t aS) \n) = 
V m J 

N J dwd9 (u a up - ^-5 a ^j C{n,w,9) \u{w,9)) + 0{l/N). 



(20) 



For example, for a diagonal element T£ (with a fixed a < m — 1), this equation gives 

(n a /N - l/m) \n) = J dwd9 (w a - l/m) C(n, w, 9) \u(w, 9)} + 0(1/N) (21) 
which in turn implies 

(w a -n a /N)C(n,w,9) w (22) 

suggesting that the function C(n,w,9) peaks sharply near w a = n a /N. 

Looking at the action of off-diagonal elements of the generators further constrains C (n, w, 9) to 
be of the form 

C{n, w, 9) w e l§ - H h(w - n/N). (23) 

Here, h(w — n/N) describes a function narrowly peaked at the origin and which satisfies the con- 
dition, 

J dw\h(w-n/N)f=(-) , (24) 

to be consistent with the normalization of the harmonic oscillator states. The detailed form of the 
function h will not matter. Note that the vectors n, w, and 9 are all (to — 1) dimensional. 
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Having determined the form of the function C(n,w,9) in equation (23), we can compute the 
matrix element between two harmonic oscillator energy eigenstates of any finite degree polynomial 
of the generators, F(T) using 



(ni, n 2 , ■ ■ ■ |F(T)|ni + m u n 2 + m 2 , . . .) 

(25) 

where the n a — > oo with the m a and w a = n a /N fixed. Roughly speaking, the magnitudes, 
w a , record the gross pattern of the n Q 's, but because the n a 's get divided by N to give w a , this 
information is not precise enough enough to distinguish n a from n a + 1 or n a + 2, etc. Those 
distinctions are produced by the 9 a dependence, which records the fine details of the states in a 
kind of holographic fashion. 



3 Low spin states 

The states that interest us for the application to large N c baryons are the low spin states. For the 
remainder of this letter, we restrict to an SU(4) spin-flavor symmetry appropriate for baryons built 
only from u and d quarks. As described in and 0, since the baryons always have exactly the 
same antisymmetric behavior in color, we can equivalently examine a baryonic state as though built 
from bosonic, colorless quarks which have the have SU(m) group theoretical properties. 
We specify our representation for the quarks as follows: 

a\ creates a u quark with spin up, 
a\ creates a u quark with spin down, 

(26) 

c4 creates a d quark with spin up, 

a\ creates a d quark with spin down. 

This assignment has the structure of a tensor product of the two dimensional flavor space and the 
two dimensional spin space, and we can take the spin and isospin generators to be half the usual 
Pauli matrices in the appropriate space. 
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In a low-spin large N baryon, there are a few valence quarks, but most of the quarks are combined 
in spin-0 pairs. If there are only two quark flavors, these spin-0 pairs also have isospin (this is 
what makes the two-flavor case so much simpler that three flavors). Thus up to a normalization, a 
low spin baryon corresponds to the following combination of bosonic creation operators: 

\N, k a ) oc z t(*-fc)/2 t*i a f2 a tk s a fk 4 |Q) (27) 

where the number of valence quarks, k = J2 a k a , is much smaller than N and finite as iV — > oo, 
and the operator 

= a\a\ — a 2 a 3 (28) 

creates a spin and isospin zero pair. We would like to write \N, k a ) as a sum of harmonic oscillator 
eigenstates, since from the previous section, we learned how to express them in terms of the eigen- 
states of the SU(m) generators. We therefore expand z^ N h ^ 2 , via the binomial theorem, and use 
equation ([Tj^) to express the expansion in oscillator eigenstates: [n = N — k] 

n/2 (n/2V 
p ]\{n/2-])\ 

■ y/(n/2-j + h)\ (j + k 2 )\ (j + h)\ (n/2-j + h)\ ^ 

■ \n/2-j + fei, j + k 2 , j + h, n/2-j + fc 4 ). 

In a typical term, both n/2 — j and j are some finite fraction of N — implicitly much larger than 
any of the k a 's for low spin baryons — so to leading order this expansion simplifies to 

n/2 

\N,k a ) oc (n/2)! ]T (n/2 - j) (fcl+fc4)/2 J (fe2+fc3)/2 

i=° (30) 

■ \n/2-j + k%, j + k 2 , j + fa, n/2- j + fa). 

To obtain the matrix elements of the SU(m) generators between low spin baryons, we can 
simplify the calculation by directly evaluating such a matrix element and then determining the 
normalization. The matrix element for a general polynomial F(T) of the group generators between 
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two baryons specified by the valance numbers k' a and k a is 
(N,k' a \F(T)\N,k a ) oc 

n/2 

(-l) j ~ e (n/2 - j)( fe i+ fc 4)/2 j( k ^+ k -i)l 2 {n/2 - £)(K+K)/ 2 (/(K+K)/ 2 
j,i=o ( 31 ) 

■ (2 £ + k^£ + , £ + k^ , ^ £ + k^ | 

•F(T)|f - j + k u j + k 2 ,j + k 3 ,^-j + h) 

In terms of the coordinates of the generators (using equations (|TU|) and and neglecting the 
0(1/ N) terms, this matrix element becomes 

(N,k' a \F(T)\N,k a ) « 

n/2 

dwdd jh (_i)i^^0-^)(^+^-ex) e ^(^') e ^(fc'-fc)/2 
j,e=o 

.{n/2 - j)(fci+*4)/2 j(k2+k 3 )/2( n /2 _ g(k' 2 +k> 3 )/2 (32) 

- 1/2 + £/iV, w 2 - £/JV, u> 3 - £/N)* 
■F(T(w,9)) h{ Wl - 1/2 + j/N, w 2 -j/N, w 3 -j/N). 

In the limit N — > oo, the sum will become an integral. Since the function h sharply peaks at the 
origin, it imposes that j/N ~ £/N, which we will call x/2 = j/N, as well as that 

Wi pa w 4 pa (1 — x)/2 w 2 Paw 3 pax/2. (33) 

The two sums over j and £ can then be reorganized into an integral over x (from to 1) and a sum 
over j — £ which becomes 

(34) 

pa 6(02 + 93-9!-%). 
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Therefore, in the large N limit, the matrix element becomes 



(N,V a \F(T)\N,k a )=A k A k , f d9 1 d9 2 d9 3 C dx 



Jo 



(35) 




+k' 4 )/2 . (k 2 +k 3 +k' 2 +k' 3 )/2 



■ F(T(w,6j) 



with 



w\ — = (1 — x)/2 , u>2 = w 3 = x/2. 



(36) 



The remaining normalization factor, Ak, is fixed by evaluating the integral for F(T) = 1 to be 



An important feature of the integral for the matrix element in equation ( |35|) is how it depends 

on the external baryons and in particular on the "valance numbers" k a and k' a . This dependence 
occurs in both the 9 terms and in the (1 — x Y kl+ki+k '^ +k i^ 2 ■ x ( k 2+k 3 +k' 2 +k^)/2 f ac ^ or w hich arose 
from the normalization of the harmonic oscillator states. Because of Bose-Einstein statistics, the 
large iV states depend strongly on the few indices that are not combined into spin and isospin zero 
combinations. This is the essential physics of the order N matrix elements of Tg between low spin 
states. 

Notice that upon setting 62 + #3 — Q\ = n, as in (|34D , the linear combinations corresponding to 
the pure spin or the pure isospin generators vanish identically. This observation amounts to the 
familiar statement that the matrix elements of pure spin or isospin operators between low spin, 
large N baryons are only of order 0(1), since the integral captures only the leading behavior in N. 

3.1 Some Sample Calculations. 

The integral formula in equation (|35D provides a simple means for extracting the leading N features 
for the low spin baryon matrix elements. We will show this simplicity with three examples. To 
begin, let us calculate the matrix element, 




1 (l + k)\ 



(37) 




(p,l/2\a 3 T 3 \p,l/2). 
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(38) 



The spin up proton corresponds to k% = 1 with k 2 = k% = k 4 = while the operator 0-3X3 corresponds 
to the following linear combination 

rpll _|_ /Ji44 _ rp22 _ rp33 (39) 

which becomes (1 — 2x), with no theta dependence. Thus 

(p, 1/2 | a 3 T3 I p, 1/2) = — — / dx (1 - x) (1 - 2x) 







(40) 



■ | 5(^ 2 + #3 " ^1 " 7T) 

= iV/3. 

This result is correct, for the exact answer is 

(JV + 2)/3. (41) 

In fact, this factor is general — the exact large N result for any matrix element between spin 1/2 
states is simply the N — 1 result multiplied by (iV + 2)/3. For N = 3, this gives the famous factor of 
5/3 for the renormalization of the axial vector current. Equation fl4~I|) (and ( fl6|) and (|49| ) below) can 
be easily obtained by explicitly constucting the low spin states out of colorless commuting "quark" 
creation operators, using (|27|) and (p8|) . One can then derive recursion relations relating the matrix 
elements for for different values of N, and solve them to obtain these exact results. But if we only 
need the leading contributions, the integral formula (|35D captures them all in a much simpler way. 
The spin +3/2 A ++ state corresponds to k% = 3, and the matrix o\T\ corresponds to 

rplA _j_ rp23 _j_ rp32 _j_ ^n41 (42) 

In the matrix element 

(A, 3/2 |ain | P, 1/2). (43) 

The phase from (^) is e~ ldl , so only T 14 (which has an e 1 ® 1 dependence) can contribute — the others 
are eliminated by the 9j integrations. Then the result is 



N / 9' / 4' r 1 
(A, 3/2|ainb, 1/2) = — ^ m ^ I dx (1 - x) 2 (l - x)/2 

■ J d6 1 d6 2 d6 3 5(6 2 + 6 3 -6 1 -7i) ( 44 ) 
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The exact result is 

y/(N-l)(N + 5)/8. (45) 

As with (f|l|), this result can be easily generalized. The matrix element between any spin 1/2 state 
and any spin 3/2 state is the N = 3 value multipled by 

y/{N-l){N + 5)/A . (46) 

One more, for good measure — let us calculate 

N 4' r 1 

(A, 3/2 | a 3 r 3 | A, 3/2) = — — J^dx{\- xf{\ - 2s) 

• J deide 2 de 3 s(e 2 + e- i -e l -Ti) ( 47 ) 

= 3iV/5. 

The exact result is 

3x(iV + 2)/5. (48) 
As before, the matrix element between any two spin 3/2 states is the N = 3 value multipled by 

(N + 2)/5. (49) 

4 Conclusions 

We believe that the integral formula, (|35|) , in addition to providing a simple calculational tool, 
yields some insight into the nature of the large N enhancement of matrix elements. The basic 
physics is Bose-Einstein statistics. The low spin large N states contain a large number of spin and 
isospin zero pairs in addition to the "valence" quarks. In these pairs, because of Bose-Einstein 
statistics, the creation operators that duplicate those of the valence quarks dominate over those 
that do not appear in the valence sector. It is this asymmetry that produces the large N matrix 
element enhancement. 
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